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Abstract 

We incorporate all existing solar neutrino flux measurements and take solar 
model flux uncertainties into account in deriving global fits to parameter 
space for the MSW and vacuum solutions of the solar neutrino problem. 



^Also Department of Astronomy. Research supported in part by the DOE, the NSF 
and the TNRLC. Bitnet: KRAUSS@YALEHEP 
2 Supported in part by the DOE. 



The solar neutrino problem continues unabated after over 20 years and 
four experiments. In spite of the relatively high flux recently reported by the 
GALLEX U collaboration, a neutrino based solution of the solar neutrino 
problem still remains strongly favored 0, at least as long as the Homestake 
CI [[I result remains unchanged. (Note that the 20 year Homestake weighted 
average used here and in most analyses is lower than an unweighted average, 
or an average over only the last 3 year's data. If a higher Homestake value 
is used, this can affect the comparative range of various neutrino-based so- 
lutions.) In this paper, we investigate what range of neutrino masses and 
mixing angles (for 2 species mixing) remains consistent with a global fit of 
all the data, including the updated 2 year SAGE average of 58^24 ± 14 SNU 
[f|]. We stress the following: (1) The MSW range is not the only range of 
masses and mixing angles which is consistent with the observations. Indeed, 
while the original MSW parameter space allowed by the Homestake data was 
far larger than for vacuum oscillations, this difference is now much less sig- 
nificant, at least for the canonical small angle MSW region; (2) There is no 
justification for incorporating the new GALLEX result and not the SAGE 
result in model fits; (3) Solar model uncertainties noticeably increase the al- 
lowed range of neutrino parameter space which is consistent with the existing 
observations. 

Neutrino based solutions of the solar neutrino problem fall into 2 major 
classes, related to suppression of the neutrino signal at the earth due to 
either oscillations between neutrino flavors or helicity states. Both involve at 
least one non-zero neutrino mass. The latter type involves potentially non- 
trivial time dependence of the solar neutrino signal over the solar cycle, but 
requires extremely large neutrino magnetic moments to remain viable, and 
has been recently investigated elsewhere 0. The former falls into two sub- 
categories: resonant neutrino oscillations inside the sun due to a changing 
electron density in the solar core (MSW oscillations ||), and oscillations in 
vacuo between the earth and sun (vacuum oscillations). The latter of these 



involves masses which are 2-3 orders of magnitude smaller, and consequently 
vacuum oscillation lengths which are up to 6 orders of magnitude larger than 
those appropriate for MSW matter oscillations. 

Both MSW and vacuum oscillations result in an energy dependent sup- 
pression of the neutrino signal. The latter has a clear oscillation in energy 
[0]. In principle, this could be useful to distinguish between these possibil- 
ities once high statistics measurements of the solar neutrino spectrum are 
available. We find that the Kamiokande || |§ spectral measurements do not 
presently constrain regions of MSW parameter space beyond those which 
can be ruled out out by global fits to overall flux measurements. Vacuum 
oscillations also allow a possible seasonal variation in the observed signal, 
as the earth-sun distance changes, and different oscillation lengths are sam- 
pled. Time sequenced data is available for each of the experiments. However, 
there is a great deal of scatter in the data, and it is not clear at this point to 
what extent any systematic time variation can be extracted from the signal 
( except that a systematic day-night variation which would occur for certain 
of the MSW parameter space is not observed fTOH). Thus, in this paper, we 
merely compare the predicted time-averaged total signals in the Homestake, 
Kamiokande, SAGE and GALLEX experiments with the observed signals 
in order to derive, by statistical methods, the region of parameter space 
which remains viable, after solar model uncertainties have been taken into 
account. Additional tests based on spectral information, or time dependence 
can serve to further reduce this region, but probably only marginally so at 
present. Several recent analyses have examined these latter issues in more 



detail, especially for vacuum oscillations JTT], [12 . 

The strategy for calculating the neutrino survival probability at the earth 
is different in the case of vacuum oscillations and MSW oscillations. In both 
cases, we consider here only possible oscillations to an active neutrino, as 
this remains at present the most likely possibility now that evidence for a 17 
keV neutrino has diminished Hl3| . For the case of vacuum oscillations H7I the 



details of the SSM production regions and N e density profile are unimportant, 
the v e survival probability P vac at a distance L from the center of the sun is 
simply 0, 



Pvac^e -+v e )=l- sin 2 29 sin 2 f ^^ J (1) 

We averaged P vac over a varying Earth-Sun distance to take into account (in 
a simple minded way) the motion of the Earth around the Sun during each 
measurement. The average was taken (linearly) from L(l — e) to L(l + e) 
with e equal to the Earth's orbital eccentricity ||. 

For the MSW model |6|, [L4| there exist good analytic approximations to 
the v e survival probability under the assumption that the electron number 
density (N e ) decreases exponentially N e ~ exp(— r/r ) near the resonance. 



We used the approximations and best fit values of tq given in |L5| who break 
the (Am 2 , sin 2 29) into several regions, giving approximations good to a few 
percent in each region. A natural distinction is between adiabatic and non- 
adiabatic transitions. If we define: 

/Am 2 \ /sin 2 29 \ 

then the non-adiabatic region is given by the condition 4no < 1 while in 
the adiabatic region 4no ^> 1 (a reasonable choice for the transition value 
is 4n = 4) [JO]]. We then subdivide the non-adiabatic region into three 
parts, depending on the size of iVJ es = Am 2 cos29/2\/2GfE compared to 
N (i) = N (o)^ + tan2 9)- 1 and N® = N^°\l - tan2^)- 1 , with JV(°) the 
electron density where the v is produced. Writing 

Am 2 
X = 27rr °^2E~ ^ 

y = 2vr no(l - tan 2 9) (4) 

cos2^ m = {l-T])/yJ{l-ri) 2 + t&n 2 29 V = N^/N r e es (5) 



the analytic expressions for the v e survival probability we used were 

N res < ^(1). 



P{y e -> u e ) -- 

N (l) < N res < N {2). 
N res > N (2). 




1 



2 1 + e- 



cos 26* m cos 29, (6) 



P(z/ e -»• v e ) = g I 1 + exp(-7rn )] . (7) 



P(z/ e -»• u e ) = - + - cos 29 m cos 20. (8) 

For neutrinos in the adiabatic region we used equation (||), which matches 
smoothly onto equation ([S]) in the limit 4n 3> 1. In addition we found that 
for the higher mass gaps including the possibility of double resonances |L6|, [L7| 
noticeably altered the results. We assumed that for neutrinos produced at a 
radius r pTOl i a fraction 
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2 . i - yi - (^ es /vod) 2 J (9) 

of them passed through two resonances. This modifies the survival proba- 
bility by Pjump —* 2Pj ump (l — Pjump) where Pj U m P is the coefficient of the 
cos 29 m cos 29 term in (|)-|) . 

Finally to compute the survival probability in the SSM we averaged (||-|) 
over the v production regions in the sun |L4] . 



For each model and each (Am 2 , sin 2 29) the survival probabilities, along 
with the SSM u e spectra |14|] were used to compute u e and v^ spectra at the 
Earth. These spectra were then convolved with simple detector efficiencies 
and cross sections 0, |], *° obtain predicted rates. In addition we gener- 



ated 100,000 "flux variations" using a Monte-Carlo procedure (as outlined in 



[0], based on earlier solar model Monte-Carlos |TJ], [T6|) to estimate the SSM 
errors in the predictions. 

The results of this stage then are a set of predicted average event rates in 
the detectors (Homestake, Kamiokande, SAGE and GALLEX) with errors 
coming from the SSM uncertainties. These predicted rates (pW ± <jW) were 
then compared with the experimental rates (e^ ± a^') by calculating the 
ratio 

(combining errors in quadrature as usual) and computing x 2 f° r a fit of the 
ratios to unity. The averaged rates [JJ, |], [3], |[ used for the fits are shown in 
table 1. 

Finally, in order to estimate the effect of solar model uncertainties and 
our fitting procedure, we compared the statistical fits described above, af- 
ter incorporating the results of our Monte Carlo over flux uncertainties with 
those obtained before these uncertainties were incorporated. Note that the 
fits for the data without standard solar model uncertainties were done using a 
standard \ 2 fitting procedure. Thus the incorporation of solar model uncer- 
tainties also changes the fitting procedure. Our results are then displayed in 
figures 1-3. As expected, the inclusion of the GALLEX result is to reduce the 
allowed MSW phase space compared to that allowed for SAGE alone. The 
inclusion of solar model uncertainties on the other hand increases the region 
allowed at the 90% confidence level for both MSW and vacuum oscillations 



(this result should be contrasted to that obtained in fTSfl). 

As indicated above, spectral measurements and a search for possible day- 
night variations can in principle further constrain the allowed regions dis- 
played here. Present spectral constraints from Kamiokande || do not yet 
provide further constraints on the MSW parameter space, but constraints 
from day-night variations do. We display in the figures those regions of the 
otherwise allowed MSW parameter space which are ruled out by these con- 



6 



straints []TIJ . 

To conclude, the progress on the resolution of the solar neutrino problem 
has been uneven. While the data does favor a neutrino based solution, un- 
fortunately at present the results of the Gallium experiments are sufficiently 
inconclusive so as to delay a definitive statement in this regard. Neverthe- 
less, when all the existing time-averaged data is incorporated, along with 
solar model uncertainties, a wide range of possible neutrino masses and mix- 
ing angles, involving either matter or vacuum neutrino oscillations, remains 
viable. We may have to await a combination of high statistics experiments 
and experiments with enhanced spectral sensitivity, such as SNO ]I|J, super- 
Kamiokande [p0|1 , and perhaps Borexino pl |, before a specific solution can 
be unambiguously determined. 

We thank S. Petcov and D. Kennedy for useful discussions, and for de- 
scribing various details of their recent work on this subject to us. L.K. and 
E.G. also thank the Aspen Center for Physics for hospitality during the course 
of this work. 



Experiment 


Rate 


Homestake 


0.28 ±0.04 


Kamiokande 


0.49 ±0.08 


SAGE 


0.44 ±0.21 


GALLEX 


0.63 ±0.16 



Table 1: The experimental rates, normalized to the standard solar model 
predictions, used in the fits. 



Figure Captions 



Figure 1: Region of mass-mixing angle space for neutrino oscillation so- 
lutions allowed at the 90% confidence level for the combined Homestake- 
Kamiokande-SAGE data, (a) without, and (b) including solar model uncer- 
tainties. Within this region the portion excluded by the day-night effect in 
Kamiokande has been removed. 

Figure 2: Region of mass-mixing angle space for neutrino oscillation so- 
lutions allowed at the 90% confidence level for the combined Homestake- 
Kamiokande-GALLEX data, (a) without, and (b) including solar model un- 
certainties. Within this region the portion excluded by the day-night effect 
in Kamiokande has been removed. 

Figure 3: Region of mass-mixing angle space for neutrino oscillation so- 
lutions allowed at the 90% confidence level for the combined Homestake- 
Kamiokande-SAGE-GALLEX data, (a) without, and (b) including solar model 
uncertainties. Within this region the portion excluded by the day-night effect 
in Kamiokande has been removed. 
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